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SUMMARY 


This report describes the analytic procedures used to prepare 
for wind tunnel flutter testing of the DAST ARW-1 -cantilevered wing 
model. The report contains a description of the methods used to: 

1) Obtain an analytical realistic simulation of the structural, 
inertial and aerodynamic properties of the wing. 

2) Obtain an analytical representation of the electro-servo 
actuator used for flutter suppression. 

3) Obtain an open loop prediction of the flutter speed at a fixed 
Mach number (.897] vs. dynamic pressure. 

4) Develop, a procedure for calibrating the flutter model at 
different dynamic pressures and to check out the feed back flutter 
suppression law. 

5) Obtain an analytical representation of a reduced state 
approximation to the full open loop simulation to be used in control 
system computer laboratory at the LaRC. 
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INTRODUCTION 


Real time, feedback control for flutter suppression is under 
serious study and consideration for aircraft (Refs. 1,2, 3, 4). The 
modal method (Refs. 5,6) is well suited for application in flutter 
suppression since the onset of flutter may be adequately described 
as a linear system instability (Ref. 7). Previous analyses of the 
problem have employed generalized coordinates, based on zero airspeed 
vibration modes or other fixed wing deformation shapes, from which 
generalized aerodynamic forces have been computed (Refs. 1,2,4), The 
contribution of this report is that physical coordinates of bending 
and torsion of the wing structure are directly employed, and that con- 
stant influence coefficient matrices are used to describe the structural, 
inertial and aerodynamic forces over a wide range of Mach numbers and 
airspeed. The aerodynamic influence coefficients were obtained through 
a modification of the SOUSSA digital program (Refs. 8,9) generated with 
the assistance of Prof. L. Morino of Boston University, Dr. E. C. Yates 
and H. Cunningham of NASA/LaRC. In contrast, the aerodynamic coefficients 
used in Ref. 1,12 were obtained using a doublet lattice method (Ref. 10). 
The method of Pade*^ approximants (Refs. 7,11,12,13) was applied to derive 
the aerodynamic influence coefficients in the real time domain. The 
structural influence coefficients were obtained through the use of the 
SPAR computer program at NASA/LaRC. The structural and geometric data 
of the DAST ARW-1 Wing, at 111 grid points, was supplied by Mr. R. Doggett 
of NASA/LaRC. 

An application of the use of physical coordinates in the real time 
domain for flutter suppression is contained in Ref. 14. The major improve- 
ment contained in this study is the application of a better state of the 
art of the unsteady aerodynamic model (Refs. 8,9). 

An earlier study (Ref. 15) for flutter suppression of the DAST 
ARW-1 wing was carried out for NASA LaRC. In order to prepare the use 
of the wing model in a flutter suppression wind tunnel test, it was neces- 
sary to develop a model of the control system actuator and integrate the 
effect of the hydraulic actuator constraint and the electrical valve 

This study was funded by NASA Langley Research Center under Contract 
NASl-15593, 
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response of the feed back control loop into the analysis. This report 
contains the necessary modifications to the analysis of Ref. 15. 

The full state simulation, contained in this report consists of 
forty-nine (49) degrees of freedom, including the actuator and the control 
surface. In order to obtain a more practical representation for use in 
control system lab work, a reduced state model is developed using the 
method of residualization. 

Finally a frequency response technique is developed, using the out- 
put of an accelerometer and a measurement of control surface rotation, 
driven by an oscillating electric current input at the hydraulic actua- 
tor control valve, as a device for checking the actuator-wing simulation 
model and to calibrate the flutter suppression feed-back law. 

An additional analytic effort directed toward reducing the number 
of degrees of freedom and the development of a control law for real time 
flutter suppression is contained in Ref. 16. 
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I . Description of the Wing Flutter Model 


The DAST ARW-1 wing was designed at the Langley Research Center as 
a swept back, cantilevered wind tunnel flutter model of a prototype, 
remotely-piloted, drone aircraft used to study active control concepts 
including flutter suppression. The data used in this report for deter- 
mining the inertial and structural characteristics of the wing, as well 
as the results of vibration and wind tunnel flutter tests, were furnished 
by Mr. R. Doggett of NASA/LaRC. 

The geometric planform and dimensions of the wing are shown in 
Fig. 1. The leading edge has a sweep back angle of 44.32°. The wing 
has a taper ratio of .392 and an aspect ratio of 6.4. The airfoil is a 
NACA 65A10 section. The main structural beam is a single tapered aluminum 
bar construction with a cruciform cross section (see Fig. 2). The dimen- 
sions of the spar cross section at various locations along the length are 
shown in Table 1. The measured stiffness distribution is shown in Fig. 3 
in terms of the bending and torsional stiffness, El and GJ curves. 

The wing is divided into eight pod sections by means of seven ribs 
oriented in the stream direction (see Fig. 4) . Each section contains con- 
centrated masses rigidly connected to the main beam to provide realistic 
mass offsets with respect to the local elastic axis. A more thorough des- 
cription of the wing model is contained in Ref. 15. 

A control surface is provided along the trailing edge, equipped with 
an electro-hydraulic servo-actuator. The surface hinge line is located 
at 80% of the local streamwise chord. The reaction torques of the actuator 
are constrained by a link to the main structural beam in the control surface 
pod section. This report contains the modifications to the analysis re- 
quired to simulate the hydraulic and structural constraints. 

A Cartesian coordinate system is used in the analysis. The origin of 
the system is at the intersection of the wing root chord and the wing 
leading edge. The x-axis is positive forward in the streamwise direction. 
The z-axis is positive down, and the y-axis forms a right hand system. 

(See Fig. 4). The structural axis of the beam defines the y* coordinate, 
rotated at an angle of 40.7° with respect to the y-axis. The origin of 
the x', y', z' system is at (-.372364, 0., 0.). 
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Thus 



where 0 = 40.7°. 


The dynamical coordinates are located along the y' axis (see Fig. 5). 
These consist of seven vertical deflections, h (y'^) (i = 1,7), seven 
rotations in the air stream direction, a(y'.) (i = 1,7), a rotation, a , 
about the control surface hinge line representing the actuator housing 
constrained by structural link to the wing beam, and a rotation, 6 , 
about the hinge line representing the control surface deflection. The 
dynamical coordinates form a 16 x 1 vector, w, defined by 


r 


h 

h 


1 

2 





C2) 


«7 


16x1 


The equations of motion will be generated in temis of the forces and 
moments affecting these 16 degrees-of-freedom. The center line, or root 
section, is constrained to maintain zero deflection in the vertical direc- 
tion and in streamwise rotation. Thus h (0) = a(0) = 0. 
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II. Structural Influence Coefficients 


The structural influence coefficients were computed utilizing the 
SPAR computer program at the LaRC computer facility. The SPAR program 
requires that the wing be decomposed into a series of grid points. At 
each point, six degrees-of-freedom are permitted. There are three 
translations along, and three rotations about each of the three axes. 

A finite element method is employed to compute the linear relationships 
between the deformations of the grid points with respect to one another 
and the resulting forces and moments resisting these deformations. For 
the DAST wing, a grid decomposition of 111 points was used, resulting in 666 
degrees-of-freedom. (See Figure 6). Of these, points 103 through 109 
correspond to the concentrated masses which undergo rigid motions without 
relative structural deformations. Point 111 is the control surface linkage 
constraint for rotation about the hinge line. 

To produce the structural influence coefficients for the 16 degrees- 
of-freedom defined by the dynamical coordinate vector, w, we constrained 
7 points in 2 degrees-of-freedom (vertical deflection and rotation in the 
flight direction) . The 15th degree-of-freedom is constrained to be a 
pure rotation about the control surface hinge line. The 16th degree-of- 
freedom, defining the control surface rotation, is constrained by the 
compression of the actuator hydraulic fluid and is computed independently 
of the beam influence coefficients. 

All other degrees-of-freedom in the original SPAR deformation are 
left unconstrained except for the root section grid points (97, 100 and 
102) which are constrained in three directions. (See Figure 7). The 
SPAR program solves a static equilibrium problem for which one coordinate 
of w is set equal to unity, and the other twenty-three are constrained 
to be zero. The forces and moments at the 24 locations are computed by 
solving a set of 24 equations of equilibrium. Thus, we have 
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( 3 ) 


where the £th element of w, w^ , = 1 and all other elements, including 
h(0) and a(0), are set equal to zero. The forces, ^ , and the moments, 

Mq ^ , are the reactions at the root section required to hold the root 
section undeformed. 

The elements of the £,th column of the influence coefficient matrix 
(K) are given by 

The deflections at the remaining 642 grid points are left uncon- 
strained. The beam matrix of influence coefficients is given in the first 
15 rows and coliomns of Table 2. The units are in Newtons per meter of 
deflection. Newtons per radian. Newton-meters per meter of deflection and 
Newton-meters per radian arranged as follows: 



Newtons/meter (7x7) 


^^Newton-meter /meter 
(8x7) 


j Newtons/rad (7x8) 

j Newton -meter /rad 
(8x8) 


( 4 ) 


To obtain the elements of the 16th row and column, we note that 
the relative rotation of actuator vane with respect to the actuator 
housing is resisted by compression of the hydraulic fluid. The fluid 
bulk modulus, 6 , is related to the change in pressure and the percentage 
change in volume by the equation 


_A 

AVo/Vo 


(5a) 


The resisting torque is related to the pressure difference across the 
actuator vane, by the equation 

2 2 
£(r/ - r/) 

T = (P^ - P 2 ) 2 
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where 


£ = length of actuator vane 

r^ = outer radius of the actuator vane 

r^ = inner radius of the actuator vane 

V = volume of the fluid 
o 

AV^ = change in fluid volume due to compression 

Since the pressure difference is twice the compression pressure 
change, Ap, we have for the actuator torque, T, the expression. 


T = 2 Ap 


where 



C5c) 


(5d) 


The change in voltime, in the actuator chamber, AV^ , is related to 
the relative rotation of the vane with respect to the housing by the 
equation 


AV = A ( 6 - e„) C, 
o ^ A 


(5e) 


Thus, we have 


2C^ B 


A(6- 


(5f) 


and 


^15,16 ^16,15 '^16,16 


V 


C5g) 


The remaining terms of ^ 26 

For the proposed actuator we have 


= 0. for i = 1 through 14. 


k. . = 19.887816 Newton meters/rad 

16,16 
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Table 2 contains the actuator housing/structural beam constraint. 


^15,15 


.81689104(10)^ 


Newton meters /rad 


(5h) 
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III. Mass Data 


The SPAR program carries out a vibration analysis of the DAST wing 
by solving an eigenvalue problem utilizing the input mass and stiffness 
data at the 111 grid points. Of the 666 degrees-of-freedom, 18 at the 
root section (grid points 97,100 and 102) are constrained to be fixed, 
and 42 correspond to the 7 rigidly attached concentrated masses (103 
through 109) . For the purposes of this study, the niomber of degrees-of- 
freedom has been reduced to 16. In order to produce a simulation in which 
the significant vibration modes are well represented, it is necessary to 
compute a set of lumped masses at the c.g.'s of the seven sections, which, 
together with the rigid concentrated masses and the 16x16 influence coef- 
ficient matrix, will reproduce the significant low-order vibration modes. 

To accomplish this, the distributed beam and plate masses have been 
summed in 7 sections and are listed in Table 3. The concentrated masses 
are listed in Table 4. 

The combined masses and mass moments of inertia (in the mks system) 
of Table 3 and 4 form the first 14 diagonal elements of the mass matrix 
shown in Table 5. The mass moment of inertia of the actuator housing is 
given in Table 5 as the 15th diagonal element. To obtain the mass moment 
of inertia of the 16 degree-of-freedom we must add to the moment of inertia 
of the control surface, m^, the equivalent inertia of the trapped hydraulic 
fluid, IgQ, The sum is listed in Table 5 as the m^^ element. 
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IV. Frequencies and Modes 


To ensure that the SPAR representation of the DAST cantilevered 
wing is a valid simulation of the wind tunnel model , a comparison of 
the first four computed and measured vibration modes was carried out 
using the full 111 grid-point model. Figure 8 contains the results of 
the vibration test furnished by R. Doggett of NASA/LaRC. Figures (9a), 

(9b), (9c), and (9d) are plots of the same modes with the SPAR program 
using all of the 111 grid points. The comparison is seen to be good. 

Figures 10 and 11 contain the next two highest modes obtained by the SPAR 
program for which no vibration data is available. 

Finally, in order to test the validity of the reduced 16 degrees- 
of-freedom model, an eigenvalue analysis using the mass, M, and influence 
coefficient, K, matrices was carried out. For the purposes of the vibration 
analysis comparison, the actuator hydraulic pressure is assumed to be 
acting, however the feed back is left in the open loop mode. The results 
of the analysis are shown in Figures (12a), (12b), (12c), and (12d), for 
the first four modes. The agreement is seen to be good. The 16 degrees- 
of-freedora eigenvalues are shown in Table 6 together with the vibration 
test frequencies and those frequencies obtained with the full 666 degrees- 
of- freedom. 


V. The SOUSSA Program and the Pade^ Approximants 

The SOUSSA digital program (Refs. 8,9) computes generalized aero- 
dynamic forces for a wing of given planform executing sinusoidal oscil- 
lations in a fixed wing deflection mode shape. The generalized forces 
are computed at a given Mach number, m, for a given non-dimension fre- 
quency, k = -^, a characteristic length, b, and a non-dimensional time 
variable T = t. The generalized forces are given in terms of the 
force per unit dynamic pressure, q. Thus we have 



j Qn 



( 6 ) 


where j = and W(jo3) is the non-dimensional deflection mode shape 

oscillating at the sinusoidal frequency, w . 

In the application we seek in this study, we desire to obtain the 
aerodynamic forces in influence coefficient form. To obtain this, a 
pre-processor was developed by Prof. L. Morino to generate 15 unit impulse 
function modes, one for each of the first 14 degrees-of-freedom of our 
w vector, plus one for the control surface rotation about the hinge line. 
The actuator housing rotation, a , is assumed to produce no aerodynamic 
forces. For the purposes of convenience we retain the W vector as a 
16x1, and the corresponding generalized force vector, L/q, as a 16x1 with 
the proviso that the 15 element of the generalized force vector is known 
to be identically zero. Thus, for the ith unit impulse function mode we 
have unit deflection for the ith element and zeroes for the other 15 
elements . 


We obtain the ith column of the desired aerodynamic influence coefficient 
matrix from Eq. (6) 


LCja))^ 

. ^ J 


i 



(m,k) + j qj^ 
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The SOUSSA coordinate system is positive deflection upward, and 
a positive angle of attack is trailing edge upward. Consequently, we 
have 


hjCj(jL))/b'l 



hyCjoj)/b 

a^Ciw) 




ayCiw) 

OpjCjo)) 

6 (jo)) 
V 


16x1 


(7a) 


By choosing the characteristic length b to be unity, we have 



(7b) 


Following the method of Fade approximants, we define a l6xl unsteady 
aerodynamic force vector (per unit dynamic pressure) , x^ , in the real 
time domain, to satisfy the following matrix differential equation 


Of) - (*p) = (« 2 > " 3) 

where the matrices F , H, , H„, and H„ are of dimension 16x16 and are to 

p 1 z o 

be determined. 

The SOUSSA generalized force is positive upward, and the generalized 
moment is positive for a positive (upward) force acting aft of the 
rotation axis. Thus, we have 

{'p 

Finally, we have from the relation between time, t, and T 
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and 






d 

dT 


W(jcj) 


(JL)2 
b ^ 



fc(jo)^ 


( 10 ) 


The resulting unsteady lift, x^Cjo)), for the non-dimensional sinusoidal 
W(ju) vector is given by the Laplace transform of Eq. (8). 

(SpCjo,)) = - (j ^ k I-Fp) ' (Hj * j -f k 

( 11 ) 

and from Eq.'s (6) and (9) we have 

(.XpCj<^)} = - (Qp(m.k))+ j Qj (m»k) (iV(ja))} (12) 

It then follows that 

(j ^ k I-Fp)'^ ('h^ + j ^ k - C^k)^ H3) = QpCm,k) + j Qj(m,k) 

(13) 

To determine the desired constant matrices, independent of time and fre- 
quency, we have recourse to the Fade'" Approximants of Ref.'s 7,11,12,13. 

In what follows below, we lean heavily on the work of Edwards in Ref. 12. 

Since there are four (4) unknown matrices to be determined, we can, at 
most, satisfy only four(4) conditions. For the first condition, we 
choose to satisfy Eq. (13) at k = 0. We have 

Hi = - Fp (m,0) (14) 

For the second condition, we choose to determine from the real part 
of Eq. (13) giving 

As k increases beyond bound, we have the piston theory limit (Ref. 17) 

«3 = -T <\.iston W 
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The nonzero elements of Opj^g^Q^j* j ’ given by (i = 1,7) 


‘^i,i +7 


2 2 
C. - C. , 

1 1+1 


- ^i+i'A 


'i+7,i ^i,i+7 


^i+7’^i+7 


12L 4 ^""- I (i - 3x + 3x 2) 

m \ 4 J o o ^ 


/c ^ - C ^ 
4a / 6 7 


(1 - xp (1 + x^ - 2x^) 


/c/ - c/\ / X ^ 1 - X ^ ^ 

<113,16 = - f ( 4 )(^ - ^ * Vl'l-"l> 


^16,13 ~ ^13,15 


4 4 3 

r ^ 1-x ^ 

4a ^6 ^7 ^1 


^16,16 


1.9431 

.8764016-. 3431794 


C^= wing chord in stream direction at the start of the ith section 


X = elastic axis in % chord 
o 


(.4231) 


Xj^= hinge line axis in % of chord (.80) 
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r 


The values of the C. are 

1 


(meters) 


.8353 


.7534 

.6784 

.6087 

.5390 

.4697 

.4205 

.3432 


We choose as our third condition to match the imaginary part of Eq. (13) 
to the SOUSSA output for the flutter-reduced frequency, k^. Thus we have 


^2 - 


F Qj (m,kp y 




Finally, following the method of Jones, (Ref. 18), we choose for F^ a 

diagonal matrix to provide stable poles for the homogeneous differential 

differential equation for the x vector. We choose for F a diagonal 

p P 

matrix whose elements are given by 


F . . 

pi,i 

p,i+7, 1+7 


p,15,15 

■p,16,16 


- av/b. 


p,i,i 


i = 1,7 


arbitrary large value (1(10) ) 
- av/b„ 


where 


b. = semi chord of the ith section 
1 

b = semi chord of control surface 

0 

a = .67034 


The values of the semi chords for the DAST wing, in meters, are given 
below: 


.39708 

.35630 

.32215 


b, = .28660 

= .25244 

b. = .21759 

6 

b^ = .18622 

b = .05049 

6 

The matrices H. , H„ , H_ and F are contained in Tables 7,8,9, and 10 
1*23 p 

respectively. It should be noted that the rotation of the actuator 

housing, 0 , (w. _) produces no aerodynamic forces. 

H lb 
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• The Hydraulic Actuator and the Hydro -Electro Valve System 


Figure 13 illustrates the manner in which hydraulic fluid pressure 
applies a torque producing relative change in rotation between the vane/ 
control surface and the housing. 

The Equation of motion for the housing is given by 


m 


15 



’^IS ^15 


w 


dt 


15 


15 

E 

i=l 


l5,i 


w. - 
1 


T 


C19) 


where 


r^^ is the structural damping produced by the relative rate 
of motion between the housing and the wing beam 

k. _ • are the elements of the 15th row of the structural influence 
15,1 

coefficient matrix 

T is the torque reaction supplied by hydraulic fluid compression 
The equation of motion of the control surface is given by 


(m 


16 




w 


16 


^^16 '^16 


T + Aerodynamic terms 


( 20 ) 


where 

r^^g the damping supplied by_ fluid leakage, valve stabilization, 
etc. designed to produce .30 critical damping (See Ref. 14) 

the effective rotary inertia of the hydraulic fluid 
bQ 

The hydraulic fluid, passing through the electro-servo valve per 
unit of time, is defined as Q^. Since the quantity of fluid must be 
conserved, it is convenient to divide the flow rate into two factors. 
The first, is the flow which accommodates the volumetric time rate of 
change due to the control surface rotation. This portion is modeled as 
an incompressible fluid flow rate and is proportional to the time rate 
of change of the relative rotation between the actuator vane and the 
actuator housing. The second is modeled as that portion of the fluid 
flow volijme which is contained in the time rate of change of the com- 
pression in the actuator chamber. This compression produces fluid 
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pressure, which in turn produces a reaction torque, T , which acceler- 
ates the actuator vane and the attached control surface. We have 


“ ^^16,16 dt ^16,16 ^'^16 '^15^ 


( 21 ) 


We now define the quantity, > to he the torque deficiency, given by 


^DEF ^ ’^16,16 *^”l6 ■ '^15^ 


( 22 ) 


Eq's. (19) and (20) now may be written as 
d2 j 15 

r — ^ w. _ + r, p w. ^ = - Z k. w. + k. , (Wi^"Wir) - (l^a) 

15 ,/l 15 15 15 - 1 15,1 1 16,16 ^ 16 15 DEF 

dt dt 1—1 


'”l6 * ^Eq’ j-1 “16 * ’^le j- "16 = ’’dEF - h6,16 t“l6-"l5’ * 

dt dt 


Eq*s. (19a) and (20a) contain the derivation of the desired stiffness 
terms for the structural influence coefficient matrix given in Eq's. 

(5g) and (5h) . 

To complete the feed back loop, we now model the electro servo valve 
as a first order command system. 


.. ... T = Y fT - T ) 
dt DEF ^ '- DEF COM'' 

where Y = - 2 ^ (Ref. 1) 

and - is the command feed back control torque. 

COM 


(23) 
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VII . Equations of Motion 


The forces acting on the wing in the air speed region containing 
the flutter speed are assumed to consist of the following: 

(a) inertia 

(b) unsteady aerodynamic forces 

(c) structural restraint to wing deformation 

(d) random aerodynamic forces due to wind gusts 

(e) a stabilizing feed back torque acting on the control response 

(f) aerodynamic forces due to the control surface deflection 


To simulate these forces, assuming small deflections, we require the 
dynamic coordinate vector, w, its first and second time derivatives. 


w and w , the unsteady aerodynamic force vector , its time derivative, 

X , the three scalar random gust variables, x,,., x-,^ and w , the 
p e> > 2D g 


actuator deficiency torque, , and feed back control torque, 


To simulate the system we define a 51 degree of freedom state 
vector, z . 

Let z be defined by 



r 




w 


16x1 


w 


16x1 


^16x1 

r 

DEF 


ID 


2D J 




51x1 


(24) 


The equation of motion for this plant is given by 


^ 51x51 51x1 ^ 51x51 51x1 51x1 51x1 ® 


(25) 
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All the elements of the vector {B} are zero except that = -Y • 

Finally all the elements of the vector {C} are zero except that 
Cj.. = 1. The scalar, w , is a variable representing the 

random gust magnitude. (See Ref. 14 for details). 

To obtain a good model for the structural damping matrix, , 
we make use of the approximation that the structure provides approxi- 
mately .5% of critical damping. 

Thus, for the case of zero dynamic pressure (q = 0) , let 



K = U (o) ^^) U"^ 


(25d) 


-1 2 

where the matrix U is the matrix of eigenvectors of M K and 
is the diagonal matrix of the eigenvalues. We have as a good approxi- 
mation of D , 
w 


D = .01 (M) U ( w.) U 


(25e) 


where is a diagonal matrix of the square roots of the eigenvalues 

of Eq. (25d) . The matrix is given in Table 12. 

To facilitate the ^analysis, we first invert the (Aj) matrix and 
write Eq. (25) in the form 


d 

at 



z S = 


I 

2 3 w 

0 

0 

0 


qM 


-1 


F +qH.,M 

P 3 
0 

0 

0 


-1 


+ A. ^ {B} . + A^ ^ ^ w 

1 51x1 COM 1 51x1 g 


M"^E 

HjM-'e 

Y 

0 

0 


m"^f 

h^m"^f 

0 

0 

Sn 
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where 


-1 



0 

0 

I 

0 

0 

0 


0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

1 

0 



51x51 


(26a) 


Eq. (26) represents the simulation of the full state. For open loop 
analysis we may set the feed back control variable, to zero. To 

check the accuracy of the model we obtain the eigenvalues of the homo- 
geneous system (with w^ and T^^j^ set to zero) as a function of the 
dynamic pressure, q ; and compare the value of q at flutter with the 
flutter results obtained by wind tunnel tests. 
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VIII. Open Loop Flutter Analysis 


The results o£ the open loop flutter analysis of the DAST wing 
using the influence coefficient method was carried out by determining 
the eigenvalues of the homogeneous part of the differential matrix 
equation (26) for different values of the dynamic pressure, q, for a 
fixed Mach niomber (m = .897) and fixed airspeed (v = 136 m/sec). 

The results of the study, shown in Figure 14, are to be compared 
to a similar plot taken from Ref. 1 shown in Figure 15. In order to 
illustrate a more detailed comparison of the pertinent flutter modes, 
we have Fig. 16 which is a plot of the frequency and damping versus 
dynamic pressure for the open loop system. 

The wind tunnel results as obtained from R. Doggett were 


m = 

.897 

q = 

5.36 kPa 

V = 

136 m/sec 

OJ = 

8.0 Hz 


The comparison is shown to be good with both the wind tunnel data 
and the analytical prediction of Abel (Ref. 1). 


IX. Flutter Suppression Utilizing the Complete Model 


The transformation of the physical state vector, z, to modal 
coordinates, u, is realized through a 51x51 complex matrix, T. 

Let, 

{z} = (T) {u} (27) 

This transformation rotates the matrix ^A 2 ) into canonical form. 

Thus , 

T'^(Ai'^A2) T = (X) (28) 

where the diagonal elements of(X)are the real or complex eigenvalues of the 
matrix (A^ ^A 2 ) . 

Following the method outlined in Ref. 14 we obtain a more conven- 
ient form of the transformation, T, in real form, we use the transformation, 
Y. We generate Y from T as follows: 

If an eigenvalue is complex, then there exists a related conjugate 
eigenvalue. Thus if 



Xi = ai + j oj. 

(28a) 

then 

^i.r “i - ^ ‘"i 

(28b) 

where 

j = 

(28c) 

If the 

corresponding eigen vector of X^ and 

are {t^} , and 

where 


= {p^} + j {q^} 

(28d) 

and 


‘'ll’ 

(28e) 

A real 

equivalent of Eq. (27) is given by 



{z} = (Y) {u} 

(28f) 

where 


{yi>= 

(28g) 


and 


{yi,i) = (qp 


(28h) 


In the event that the ith characteristic root of (A^ ^^2^ 


real, then we choose 

the eigen vector so that 


II 

•H 

(Pi) 


(28i) 

The corresponding elements 

in the canonical real Jordon matrix 


( 00 ) are given by 




If the roots are a. ± oo. . 

then 



1 1 ' 



00. . 

= a. 



1,1 

1 



OJ. . 1 
1,1+1 

= 00. 
1 


(28j) 

00. . 

1+1,1 

= -00. 
1 



1+1 ,1+1 

= a. 
1 



and if the ith root 

is real then 

we have 


oo. . = 

1,1 

a. 

1 


(28k) 


The eigen vector transformation listed in Table 13 for 
the base case of Mach number = .897 and a dynamic pressure of q = S.OkPa. 
The corresponding eigen values are given in Table 14. 

The modal form for the full 51 DOF system is given by 


^ {u} = ( 0 )) {u} + (Y)'^Aj)'^ (29) 

+ (Y)'^ {C} Wg 

Examination of Table 14 indicate that the first two complex eigen 
values -have real positive parts and are unstable. Fig. -(14")' indicates that 
the open loop flutter q occurs somewhat earlier than q = 8.0 kPa at a 
dynamic pressure of ^pujYTER ~ suppress the flutter con- 

dition we design a feed back gain for T^^^^ to place the unstable pole 
at a desired location in the stable, negative real half plane. The 
method for accomplishing this is developed in Ref. 14 and described herein. 

' 25 


Let the unstable pair of coupled equations (u^ and U 2 ) of Eq. (29) 
be given by 


d 1^1 

Xl“iy 1^2 C2 


T 

COM 


(30) 


where is positive. In the above we may neglect the random gust wg 
since this external force plays no role in the unstability. We choose 
the command torque, have a fixed constant value over a fixed 

interval, AT, and equal to a linear combination of u^(tj^) and U 2 (t^), 
where t^ is the start of the interval AT. Thus, we have 



Ai 

\^21 ^ 2 y ^^ 2 *^^!^ 


(31) 


where the constants, f^^^, f^ 2 > ^ 21 * ^22 determined. It 

is plain that the feed back gain is given by 


COM 


17 (hi ^ h2 “2 'h’ ’ 


or 


COM 


(f,^ u^ (t^) + f22 «2 


^2 "21 1 ' 1 


(32) 


Thus only two of the four f . . 

11 

Under the assumption of 


are independent. 

Eq. (31), Eq. (30) may be written as 


d / ^1 \ ^l\ / f^W ^12\ / ^1 ^^1^ \ 

dF |U2(t)| Uo)^ ^ V21 ^22/ A2^^lM 


(33) 


Since the non homogeneous forcing function is constant we may integrate 
Eq. (33) and have at the end of the fixed interval 


where 




AT 





^1 

- 

e 

cos (o),AT) 






a. At 





^2 

“ 

X 

e 

sin ((jOjAT) 




^21 

= 

-^2 





CM 

CM 

= 

^1 





^11 

= 

S^2 

- aj(l-b^j) 

)/ 

r 2 

(Oj 

2- 

+ (jJj ) 

^12 

= 

(aib^2 


)/ 

r 2 

(a^ 

2a 

+ Wj ) 

^^21 

= 

-"l2 





^22 

= 

Si 






(34a) 


We now define a unique, equivalent 2x2 system, with poles in the 
prescribed location (a ± j 6 and a > 0) such that 


_d 

dt 




Uj(t) 


U2(t) 


(35) 


and choose 


® 11 ’ ® 12 ’ ® 21 ’ ®22 


®11 ^ ®12 = 20 > 0 


®11 ®22 ■ ®12 ®21 0+6 


(36a) 


and demand of the output at t = t^ + AT, that if 

Ui(ti) = u^(tj) and if U2(t^) = ^^(tp 

then at t = t^^ + AT we require 

Uj(tj+ AT) = u^(tj+ AT) 

U2(tp AT) = U2(t^+ AT) 


(36b) 
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The output of Eq. 35 is given by 


u^(t^+ AT) 
U2(tj+ AT) 


= [x^ I + X2 E] 


Uj(tp 


where 


X 


X 


1 

2 


e ^^"^(cos (6AT) - sin (SAT) ) 

aAT sin(6AT) 

® 6 


( 37 ) 


(37a) 


The two systems (Eq. (33) and Eq. (37) are equivalent in the sense that 
they periodically coincide at every AT interval provided that 


B + 

DF 

= x^ I + X2 E 

(38) 

The 8 unknown 

constants in F and E required to satisfy Eq. 

(38), 

(32) and Eq. 

(36a) are given by 


®11 

= 

K a2 - 2a - bjj - b^2/ 


®22 

= 

2a - e^j 


®12 

= 

bi2/x2 * (Xj * Xj 632 - bjj) /(V Xj) 


®21 

= 

* '*1 * ’'2 ■ hd '"'’‘2 

(39) 

^11 

= 

-r cxj t X2 - bjp / c ^ bjj * bjj) 


^21 

= 

.2 

hi 


^22 

= 

h 

* ’‘2 ®22 - - '’12 


^21 

= 

h 

= f 

£2 22 
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where 


2 x2 

K = a + 6 


^ = (^2 ^11 “ ^ 12 ^ ^ *-^2 ^12 ^ ^1 ^ 11 ^ 


(39a) 


The utilization of this form of the stabilizing feed back torque 
will in general produce a stable system with the poles of u^ and U 2 
at the desired location. The remaining poles of the system will not 
remain invariant and some shifting will occur, although the entire 
system in general, will remain stable. 

The actuator feed back law requires an estimate of the model state 
u^ and u^ at a periodic sequence of times t^, t^ + AT, t^ + 2 AT, etc. 
Since the model state is not directly observable, an estimator (or 
observer) will be required, based on some known filter method (such as 
Kalman or complementary) . Utilizing the output of some measurable 
quantities, such as an accelerometer placed at some fixed wing station, 
we may proceed to estimate the desired canonical states, Uj(t) and U 2 (t) 
at the periodic intervals. 

For the purposes of this study, which are to predict the response 
of the actuator-wing combination at various values of dynamic pressure to 
calibrate the digital feed back law, we will assume that we possess 
perfect knowledge of the 51 degrees of freedom (neglecting the effect of 
the external random wind gust , wg) . 

We now proceed to develop an analytical procedure for testing and 
calibrating the actuator-wing combination and to determine instrumentation 
for such a purpose. 



X. Testing and Calibration of the Actuator Feed Back Law 


Standard procedures for testing and calibrating wind timnel models 
are well known. Static deflection tests to obtain checks on the struc- 
tural influence coefficients, as well as vibration "shake" tests to 
obtain the normal modes and frequencies in still air are typical of these 
methods, A somewhat more sophisticated calibration and test check-out 
procedure is to use the frequency response of some measurable output of 
the wing (such as an accelerometer, or a surface rotation) driven by a 
sinusoidal input at a given amplitude and frequency. This method has the 
advantage of incorporating the entire system including the interaction of 
the structural, inertial and aerodynamic subsystems, into a single measur- 
able quantity and which has the advantage of being computable and pre- 
dictable. By recording the accelerometer and control surface rotation 
amplitude and phase with respect to the known oscillatory input signal at 
the actuator valve and comparing this with the predicted amplitude and 
phase, based on our analytical simulation of the wing/actuator combination, 
we will be testing the accuracy of our simulation and at the same time 
calibrating the feed back law. To determine the vertical acceleration of 
the wing at a point (X, Y) , let 

Yi ^ y ^ 7i+l ’ C40) 

then, the acceleration of the wing h(x,y) is given by 


h(x,y) = h.(x.,y.) + ^ ^ " 

^i.l - ^i 


(aiCx.yi 


) + 


“i+l^Vryi+1^ -a\(x,,y,) 


^i.l - >^i 


(y-ypjCx-yg^^ (y)) 

(41) 


where 


= -.372364 -.860y 


(41a) 


Thus, the value of the acceleration at any point along the wing can be 
determined from a linear combination of at most 4 coordinates of our 
physical state. 
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We choose for the simulation of the frequency response of our 
51 degree-of- freedom system the modal description of the state (Eq. (29)). 


(u(t)} = C«) (u(t)) .(« J . a 


(42) 


where, in Eq. (29), we have set 


w =0 
g 


and 


(Y)'^ (Ap-^ (b.) = (£.) 


(42a) 


(42b) 


We choose for the feed back law, the gains obtained in Eq. (32), 


'^COM £. *^^11 ^1 ^^1^ ^12 ^^2^^!^ ^ 


(32) 


or 


COM 


£. *-^21 “l *-^ 1 ^ * ^22 ’^ 2 *-^ l ^ ^ 


The scalar forcing function, a is the oscillartory input to 

the actuator valve with fixed amplitude, a, and frequency, w. (j = /^T) . 

We seek to obtain a dynamical system, equivalent to that represented 

by Eq. (42) in which we replace the piecewise constant feed back control, 

^COM’ ^ continuous function of u^(t) and ^^(t) which coincides with the 

true T„_., at the time intervals t , t + AT, t + 2 AT, etc. To accom- 
COM o’ o 0 

plish this we note that we may invert the 2x2 transition matrix of Eq. (37) 
and obtain 


Uj(t;^) 

UjCtj) 


Uj(tj^+ AT) 


/^®11 

c (-AT) I + X (-AT) 

\21 ^22jj\^2^^1* 


(43) 


It follows that in the equivalent system 
1 


COM 


£ [^11 ^12^ l^iC-AT) I + x,(-AT) 

i 


®11®12 


■® 21 ®22 



(44) 
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Our equivalent, homogeneous system becomes 




C45) 


where 

(E) = 

(P) 49 x 2 is a 49x2 matrix 

and ^^^49x49 bottom of our original real Jordon form of Eq. (42). 

The elements of (p) are given by 

a. 

Pi,l ^^11 ^2 ®lP ^ ^12 ^2 ®21^ 

£. 

Pi , 2 "" *^^11 ^2 ®12 ^12 ^2 ® 22 ^ ( 45 b ) 

(i = 3, 51) 



We may now obtain the frequency response of the equivalent modal state 

u(t), to the periodic forcing function a Since all the roots of 

the equivalent system (Eq.(45) ) are now stable, we may obtain the steady 

iw”t 

state frequency response to a sinusoidal input, a e , by 



To obtain an explicit solution for the inverse of the matrix. 


we have the following: 

Let the matrix A(NxN) be given by 



(N=51) 


A (NxN) 


(47a) 
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Let the characteristic polynomial of C AI-A) be given by 


I X-A 


N 

Z 

k=0 


,N-!k 


(47b) 


(^o = 


Let the inverse of (IX -A) be (C) . Then (C) is given by 


(C) 



1 




(47c) 


where the matrices B, are obtained by recursion 

k 

B = a I 
1 o 

(47d) 

\ ' “hc-i I * *\-i 

(k = 2,N) 

We have that every element of the inverse is a rational polynomial 


Cij(A) 


N 

Z 

k=l 

N 

Z 

k=0 


X^-^(b ) 

ij 


a, X 
k 


N-k 


(47e) 


where (b,) is the i,jth element of B, 

^ iJ 

It follows that the time derivative of the ith modal coordinate is given 
by 


_d 

dt 


u^(t) 


(jwe^'^^) 


51 

Z 

k=l 


Ci^k Cjw) 


(48) 
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we have recourse to 


To obtain accelerations ^i+1’ “i “i+1 

the relationship in Eq. (28f) 

{z(t)} = (Y) {uCt)} (49) 

so that 


h.(t) 

51 
= E 

j=l 

^i+7,j 

Uj(t) 

hi,l(t) 

51 


Uj (t) 

a^(t) 

51 
= E 

j=l 

^i+23,j 

u(t) 


51 
= E 

j=l 

^i+24,j 

u(t) 


Using Eq. (41a) we may obtain the amplitude ratio and phase for:the pre- 
dicted acceleration response to the oscillatory input. 

Another measurement which can readily be made is the control surface 
rotation, 6, which is given by the sixteenth element of the {z} vector. 
We have 

51 

<S(t) = z^^Ct) = . yi6,i ^i^^^ 


where 


u^(t) 




Cjw)^., 


Jwt 


(51a) 


Utilizing the measured and predicted frequency response for an 
accelerometer and the surface rotation for 3 or 4 values of dynamic 
pressure (below the flutter q) will serve to test the adequacy of 
the simulation and to calibrate the feed back actuator law. 


XI . Flutter Suppression Utilizing a Reduced State 


For the purposes of developing a digital discrete flutter suppression 
feed back law and to test the electrical circuitry in a control system 
laboratory it is often convenient to have a reduced simulation model of 
the full 51 degree-of-freedom model. The most favored technique for 
state reduction is the method of residualization. In this method, the 
state is divided into two categories, the slow/low frequency states that 
are most active in the flutter instability, and those states where 
natural frequency is too large to affect the flutter instability. The 
sense in which the fast states are to be neglected, or eliminated, is 
to set the time deviation of the fast states to zero and to use those 
differential equations, with zero derivatives, as a means for solving for 
the fast states as functions of the slow states. In this way, the effect 
of the fast states on the slow states is maintained in the remaining 
differential equations and the system is reduced. 

For the purposes of arranging the states according to frequency, we 
choose the real canonical form of this system simulation, Eq. (29). 

Let the characteristic roots be arranged in the order of increasing 
absolute value 


a. + 1 w. 
1-^1 


^ l^i^l ^ j 


w 


k+1' 


(i = 1,51) 


(52) 


We choose to include in the analysis all canonical states whose 

absolute characteristic amplitude is less than five(5) times the flutter 

instability characteristic magnitude. We divide the canonical state 

(u. } into two states 
1 




( V 

( fig) {Ug} 


.{£g}T 


COM 


COM 


{C } w 
A g 




(53a) 


(53b) 


/ 
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We set 


-Tt W = 0 


(54) 


and from Eq. (53b) we obtain 


'V = tV^'V^COM" 'V' 'V“g 


(55) 


The transformation (Y) , which carries {u^l into the physical state 
(Eq. (28f) ), may be decomposed into 


and 




(56a) 


(56b) 


Thus, we may represent the entire 51 physical state vector in terms 
of the reduced cononical state {u^}(Eq. 53a) and the non-homogeneous 
forcing functions and w^. 

It should be noted that the real Jordan matrix, ( fi), is a function 
of the dynamic pressure, , for a fixed Mach number, M. In carrying out 
the reduced state feed back law development, a sequence of ( (q) ) 

matrices and associated real characteristic transformation (Y(q) ) must 
be processed, so that the digital feed back law will be a function of 
q for fixed Mach number, M. 
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TABLE 1 SPAR GEOMETRY DETAILS 


DISTANCE 
ALONG 
ELASTIC AXIS 
METERS 

A 

METERS 

B 

METERS 

C 

METERS 

t- 

METERS 

0 





. 17145 

. 0331216 

. 0722122 

. 212344 

.00508 

. 393192 

. 021082 

. 059436 

. 201422 

. 00508 

.78232 

.01905 

.055372 

. 182372 

,00508 

1. 117346 

.017018 

.050038 

..17526 

.00381 

1.452118 

. 015748 

.044196 

. 150368 

. 00381 

1.787398 

. 014097 

.03937 

. 128524 

, 00381 

2. 122424 

. 013208 

. 03302 

.10287 

. 00381 

2.426208 

, 010922 

. 031242 

. 094488 

1 

. 00381 


B 



Fig. 2 Sketch of spar cross section along elastic axis 
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Figure 14. ‘ Open loop 'eigenvalues of ‘flexible modes as a - • 
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TABLE 2. MATRIX OF INFLUENCE COEFFICIENTS, 


ROW 

COL 1 

2 

COL 3 

COL 4 

COL 5 


1 

.22677542000000E+08 

-.70437274000000E+07 

•23725273000000E+07 

-.453629640000006+05 

.910847530000006+05 

2 

-.70437274000000E+07 

.66060031000000E+07 

-.43556373000000E+07 

.142402110000006+07 

-.28 59064 7000000': + 06 

3 

.23725273000000E+07 

-.43556373000000E+07 

.50189197000000E+07 

-. 2961 0165000000 E+07 

.961076810000006+06 - 
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-. 45 36 2 96 4 000000 E +06 
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. 22799027000000E+04 

-.71 5 6398 7000000 E+04 
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-.25493S03000000E+06 
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.753147080000008+03 

9 

.92075518000000E+05 
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-.358658550000006+05 

.223953280000006+05 

-.449315310000306+04 

10 

-. 33282973000000E+05 

.39630705000000E+05 

-.957801750000006+04 

-.209672530000006+05 

.1409268'^OOOOOOE+05 

ll 

.55316836000000E+0'4 

-.17384829000000F+05 

.2022>’7640000006 + 05 

.582421670000006+04 

-.21338244003000F+05 

12 

14706 400 OOOOOOE +04 

.46161998000000E +04 

-.155253920000006+05 

,?259?046000000E+05 

-.841797470000006+04 

13 

.48248703000000E+03 

-.151448160000006+04 

.509060710000006+04 

-.17144015000C00E+05 

.339447650000008+05 

14 

-.29403225000000E+02 

.92293970000000F+02 

-.310226400000006+03 

.104438190000006+04 

-.331506650000006+04 

15 

■ -. 13415067000000E+03"" 

.42108641000000E+03 

-.141539160000006+04 

,476493940000006 + 04"“ ' 

-.151343400000006+05 

16 

0, 

0* 

"07 ■ 

0." ■' . - - - 

0. 


R0l< cat 6 'COL 7 ■ COL 8 ' ’ 'C0L"9‘ • — COL'lO" 


1 16839310000000E+05 ,r2799027000000EV04 -,25<»93"86‘3T»00005E + 06 . 9207551 BOOOOOO E + 05 33282973000 oOOE*05 

7‘ . 528570090010000E + 05' 7 1 5 63'98 7000000E+0'it -.l‘tl30021000000'E + 05 - . lO'^97fl270OOOOOO E+05 ^3 9630705b0O'0b0= +05 
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TABLE 3 


COMBINED BEAM AND PLATE ELEMENT MASS PROPERTIES 


Section 

No. 

Mass 

(kgm) 

I 2 

(kgm-m ) 

2 

(kgm-m ) 

(kgm-m^) 

. % 2 
(kgm-m ) 

2 

(kgm-m 

* 

0 

7.1088 

.0262 

.0654 

.0913 

-.0487 

.0429 

1 

2.5310 

.0048 

.0343 

.0382 

.0218 

.0173 

2 

1. 8328 

.0026 

.0199 

.0220 

.0126 

.0100 

3 

1.3659 

.0017 

.0128 

.0142 

.0081 

.0064 

4 

1.1339 

.0012 

.0107 

.0116 

.0066 

.0052 

5 

.9982 

.0010 

.0095 

.0103 

.0059 

.0046 

6 

,7266 

.0002 

.0058 

.0061 

.0034 

.0026 

7 

.3341 

.0001 

.0023 

.0024 

.0014 

.0011 


16.0313 


* 

Wing root section mass properties 
** Compares with 15.8773 kgm calculated by SPAR 


TABLE 4 


CONCENTRATED MASS PROPERTIES IN X-Y COORDINATE SYSTEM 


Section 

Mass 

Pitch 

Yaw 

Roll 

** 

AX 

AY 

Pitch 

No. 


Inertia 

Inertia 

Inertia 



Inertia 


(kgm) 

2 

(kgm-m ) 

2 

(kgm-m ) 

Sc 

2 

(kgm-m ) 

(m) 

(m) 

lY(grid point) 
2 

(kgm-m ) 

1 

1.6556 

.0734 

.0801 

.0134 

-.0051 

+.0152 

.0734 

2 

1. 0342 

.0430 

.0443 

.0044 

+. 0147 

+.0053 

.0432 

3 

1. 1975* 

.0333 

.0344 

.0042 

+. 0239 

+. 0036 

.0340 

4 

1.1612* 

.0267 

.0280 

.0034 

.0201 

.0028 

.0272 

5 

1.0161 

.0177 

.0192 

.0027 

.0218 

.0041 

.0182 

6 

1.4198 

.0200 

.0206 

.0047 

.0285 

.0000 

.0211 

rj*** 

.6350 

.0101 

.0127 

.0021 

-.0521 

-.0112 

.0118 

y *** 

.5534 

.0056 

0 

.0056 

.2070 

-.0210 

.0293 

y !)<♦* 

,4627 

0 

0 

0 

0 

0 

0 


9. 1355 

I^, . , . ^ = I + m AX^ 

Y (grid point) Y 

% 

These values were altered to agree with SPAR computer model 

AX - Mass offset from grid point in X direction 
AY - Mass offset from grid point in Y direction 
*** 

These masses were added to cause flutter within the available wind tunnel dynamic pressure. 










TABLE 6 


COMPARISON OF VIBRATION MODES 



Vibration 

(Hz) 

Reduced 


Test 

SPAR 

Model 


Freq. 

666 DOF 

16 DOF 

1 

5.22 

5.279 

5.432 

2 

19.44 

18.90 

17.532 

3 

26.04 

26.01 

24.76 

4 

45.96 

44.34 

41.38 

5 

N/A 

60.19 

45.73 

6 


73.87 

60.64 

7 


84.96 

68.85 

8 


93.18 

91 -.34 

9 



107.29 

10 



125.08 

11 



135.51 

12 



147.18 

13 



203.88 

14 



296.70 

15 



420.03 

16 



894.13 


62 
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TABLE 8. 

MATRIX. 



ROW 

COL 1 

CQL 2 

COL 3 

c ol 4 

COL 5 f 

1 

1 
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13 
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14 
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15 

0 • 

‘0 . - . - 

■ 0. 

o; • ” 

■Q." ■ 

1 6 
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COL 9 • 
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15 
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0. ‘ ^ 

'o; ■■ ' " ■ 
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0. 
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TABLE 9. Hj MATRIX. 













TABLE 10 


ELEMENTS OF DIAGONAL MATRIX F 


-.22960E+03 

-.25587E+03 

-.28299E+03 

-.31811E+03 

-.36114E+03 

-.41898E+03 

-.48956E+03 

-.22960E+03 

-.25587E+03 

-.28299E+03 

-.31811E+03 

-.36114E+03 

-.41898E+03 

-,48956E+03 

-.53875E+07 

-.72230E+03 


TABLE 11. Qj^ (m,o) STEADY STATE LIFT AND MOMENT. 


COL 2 


COL 3 


COL ^ 


COL 5 



,<t85 35683000000E+00_ 
r3'5 1 V 353 230 00 0 0 EV 0 0 
1^2 2 1 1 V9 6^0 000 00 E »06 
, i 5 661 3'50 0 Ob d 0 0 E + 0 0 
;ii663032000000E-*-dd~ 
^7505614 100d0d0E-01_ 
,3915147 600 0 00 O'E -0 1 
,'617 77 233000 000E'-01' 
,‘4812741100000 0E-0^1_ 
7z 9 94 o 1 i'4 Ob o b 0 0 E-b i 
,1878912f)OOOOOOE-Ol 
;12130249000000E-OY_ 
r76la4509000000E-02 
,42734592000000E-02 


31310009000000E+^^ 
310762'93000000 E'+00_ 

2 372' 2 4 7 7 0 0 0 000 E'+ qO_ 
16 2' 3 5 2 5 9 0 Ob 0 00 E ♦ 0 q_ 
11066954000 ObOE ♦ 00 
'6066 30660 00 Oo6 E -0 1 

'3 5'2 0 9 6 9 5 0 00 0 0 OF -q^;L 
22427634000000E-0^ 

3 6'3 7 5 a 9 4 0 Oo 0 Ob E-O l 
3"1 6 5 6 5 0 1 0 Obb 0~0 F - 0 1_ 
2004R2 240 0b000E^-^'l 
1 ^6' 3'6 2 8 0 0 0 0 0 0 OF -0 1_ 
69674107000000E-02 
3832345?OOOOnO'E-02‘ 


r4O396O0TD 00000 E -0 4 


,357751160OOO0OE-O4 


.14552497000000E+00 
;2O5O9O50OOOOOOE+OO 
^^347077000000_0|_+0q' 
. 1 9976 137 0~OOOOOE+O'0j 
^1 3^2'5 8 9 3 eqO 0 0 0 0 F^+ oq_ 
^3 015 8 1 4 Oo cm 0 0 F ^ 0 1 _ 

, 3 5'995'660 000 00 0':-dl 
^ 4 2'5 6J 1 8 500 0 0 0 OE "fO 3_ 
; 1 4 3 0 7 8 8 5 0 Ob 0 00 F - 0 i 
. 2617i677bbOOOOF-01 
; 23 618 39 Ibj) 0 00 Ol '-m 1 _ 
n 42 0 1 4 4^40 ObO 0 OE-O'l 
,7 4 i 8 5'3 'i 3b"0'd0 O'bj - 0 2 
. ^89 61 3 Ibq O0 0bbF:-02" 
• 

r3T^T4o’d2 r00 90DlP“^^ 











TABLE 11 



(CONCLUDED) 


coins CoL 14 ■ COL 15 


7T4^3357^TOO^015Wyr=n5'2 r?T5V2 V9JOirOiyO O E -O 2 07 

7TD 6 V5 6 39 000 O'O'O E=OT T3'6 67 0 0 7 <ro ODWOE -0'2 ^07 

72TrO 5 62T0 0 00 0 OT^ 75'9 8 7 V8 0 7 OWOWE -02 07 

^ 5'n V9 05 200 O'O 00 E =0T TTBO 76891000000 E -01 ^07 

711532936000000 ETOO 732610500000000E-01 ^07 

717563 30000O000EV0O 787129017000000E-0I ^07 

: 8 ft 8 6 5 9 5 7 0 0 0 0 0 0 E - 0 1 1 1'2 9 8 7 6 5 9 0 0 0 00 0 F+0 0 07 

; 3 85 95460000000 E-0 3 . 71 786065000000 E-0'4 07 

763‘462952000 000E-03 ;'92 297 393000000C-04 0. 


.1"31 95 379000000 E-02 73042 1 344000000E-03 67 

.7'2405774000000E-0‘3 , 915 19052000000T-0'3~ 0. 

' -.5 42 59031000000F-02 ;T47 994 57000000 E-0'2 07 

-.15158469000000E-01' ^26361 4 89000 000 E-02 0, 

~=7992 9_4«7 3 0 00 000 E -02 -7115 8 5 3 5_2 0 0 0000 E-OT 07 

0 « 0 * 0 « 

; 39024026000000E-03 7 1 85 55 74‘5000060E=0 3 07 




N) 


TABLE 12. MMPIfJG' COEFFICIENT MATRIX D, 


W 


ROW 


COL 2 


COL 3 


COL 4 


COL 5 


9<»23e33C450974E*C2 
104317P209OB<V<<E + 02 


,36055563170l5LE*0l 


,330'3Qn60316027E + C0 


locsRaciiceaASE+co 


263256P24F1 4 74 E-01 
111 3 3137E-02 
P ?eS14732IP A 1 2E4Q0 
2B546QP2533727E+1' 0 


^41O01 3FA0fr241F-Cl 

7216515401Q545F-C2 


30150196372q2F-t 3 


.16075U92P4277E-O2 


,PP60f)3P2')l'^3F4.)[-C 2 


133S97<5F‘V6'?423E-03 


ragE!gW 


•.19431782099851E4Q2 
. 3 60 462 50709902 E 402 


- 17049Be3074519E+u2 


2356643091127PE+01 


•.310516e7170430£4QU 

.134966P620765JE-01 

.206716479652O6E-01 

•.23532259fl5e794E400 

-.409540126643fi6E-OX 


.1Q981P3P5P7123E4C0 

-.23554660323266E-01 


.83811849u4.,t01F-02 


-.9 6902917834391 E-02 


.34758068794919E-01 


.1P271507426827E-03 


- 40411051328436C-01 


.3605556317Ql^6g»01 

1704 9883 074 51 9E402 


.2e55531in09563F+02 


-.1330098‘^47753 9E+02 


, 16967617746 ~F P9e4-01 
-.3 2C23946243911E4-00 
.9489 9 859 7P« 2Q6F-01 
.590 35 76'F5 4777 »c-01 
-.187880587090378+00 


-. 113068402113298-01 
.1414 77PO<}3Q609E+00 


-.2950764P510172E-01 


- 715624125O415P.E-02 


lfc6364473423FE-01 


-.158863C.21360068-C2 


- 60235207763070E-01 


-.33099060316053E+00 

.23566430911245E+01 


-.13300985477539E+02 


.21418814072643E+02 


-.1044400C 41 4250E+02 
.11 3 63 991829943E + 01 
-.28 01 4 5 80 06608 lE-01 
•-1 15 3JP 6231 293 0 1 E-0 1 
.747'9719983 62e6£-01 


-.11 3 8e749622210£+ 00 
.49o92533j5'9373E-01 


.12919695575782E+00 


-.43640320356747E-J1 


.31157577162947E-01 


.5t52e29C652692E-02 


- 32940443848187t-01 


100892U1109063E+00 


3105168717C840E+00 


16967617746817E+01 


10444000414256E+02 


16299 22899 fc0buE +02 
79130S 5 B 262166E + U 1 
1 0280 7~239 8_C4 2 6 E +01 
7899 i'434 13 C 0 0_6 E £o 3' 
1138462512261 1E-6T 


2917 PIS 24 b E - 01 
15i/144 4420 56e5E+o0 


479454591F3341E-01 


12922653243488E+0O 


324214375S7o7bE-01 


226869lO5oeu55E-01 


175446625ot436£-Oi 



COL 6 


COL 7 


COL 8 


COL 9 


COL 10 


1 

-.263256824E1683E-1 1 

.75131453463362E-C2 

-.858514732183938+00 

.2854699Z533707E+00 

-.541901364060648-01 


2 

.134966P6204253E-01 

.2C671647965129E-C1 

-.23532259858743E+00 

-. 40954012663654 E-01 

.199618385e7208E+00 


3 


TT 1 n n Ti-HiHJa+IW 

.590357675482668-01 

-.187P805P709946E+00 

-.113066402099496-01 


4 

.11363991P2O913E+01 

-.28014580065410E-01 

.1153P623129721F-01 

.74797199837C67E-01 

-.113E8749622i;63E+00 


5 


.i028C723O8u414E+01 

-, 789914340964 17E-03 

-.113846250327046-01 

.281920286197x66-01 


6 

.120131100471998+02 

-.4729149u201940£+ul 

-.156141832768635-01 

-.362453430484476-01 

-.437392995204458-01 


7 

-. 4729149020 1949E +01 

.37457446772104E+01 

.23744566173P5PE-02 

. 5 89 60 625 8 76709 E-02 

.89697699464620E-02 


8 

-.156141P3276894F-01 

.2374456 6175321E-02 

.7ecu412i687595F+no 

-.853299376357256-01 

-.673604279453536-02 


9 

-.362453430472758-01 

.5P96062587P123E-02 

-.853299376357208-01 

.35613731413008E+00 

-.759784644105616-01 


10 

-.43739299519895E-01 

,89fi97e99464488E-C2 

-.67360427945RC3E-02 

-.75978484400637E-01 

.265405642201426+00 

1 

11 

-.20342631R129A2E-02 

.10021353165C41E-01 

.12P00587052473F-02 

-.330523775674636-02 

-.509950729998766-01 


12 

-.8235466e7374RjE-Cl 

.33555064798169E-01 

,4647340460371 28-03 

.61621401948467E-04 

-.56667103577421E-02 


13 

-.494349496371 29E-01 

-.26P30536561611E-C1 

-.199929668717408-02 

-. 52i551iE925262E-02 

-.706984515667966-02 


14 

,247250266687478-01 

-. 6678 76 38 2o6P 158-01 

,75f 86213018P1 48-02 

.17796527109582E-01 

.186847617626716-01 

1 

15 

-.447095514375938-02 

.234503193591806-01 

-. 407746563 Q3779F-04 

95.^440485316576-04 

-.150096475144226-03 


16 

.7CO23O36844029E-11 

-.4 3060300757901 E-02 

-.923041663425018-02 

-.21990114318800E-01 

-.243674082776406-wl 

1 





TABLE 12. (CONCLUDED) 


COL 13 


COL !<♦ 


COL 19 



»7216515^9gl88<?E-C2 
235546603 2?< >9nE-01 
^lA77F,993°7O9F*O0 
.<IQQ92533C.60‘t23E-t 1 
•.1501A<tAA20!:feR7E-*-00 
•.203A2631P23<.A0F-G2 
.100213S31fc5815E-Ql 
.12eC05P70S3I05E-( 2 
Is ?10 9?3775f fr?7^F-C2 
509 0 5 Q7?999e?2F-Cl 
2 <.'207^7272 2F-»C0 
-.<>2 25^ QfC. 0< P.65 3 F-C1 
-.1174 ^ AE396A45r-C 1 
.1253'o707162<-79F-01 


-.71<.96<.25P372-j 3F-C-4 
- 192«66<.615115iF-Cl 


•.81015019609527E- 
.R3PllPA90<t4367E- 
•■295U76A85095A9E- 
.129X9695575e5t/E-t 
•.<t79A54591P3516E- 
-.P235<i66P.73P22AE- 
.33555O64796053E- 
<^6<i734Ci4bl67PXE- 
61621^0207S265E- 
•.SbPb7T03575Po9E* 
-.<i225AO50066570E- 
.127394,O2210793E-i 
-. 33C29917366326E- 
- 217Q 5 Pb<te5920<iE- 
702<i52A1918493E- 


670429993P3535E- 


■.160751592B3304E- 
•.9b902917P3321bE- 
•.715b2^12?Q3’'4QP- 
-.43P4032035b6«5E- 
.129^2e53?A3453E^ 
•.A9A3<>9«t9b371fe7E- 
■.2663053^561*^139- 
•.19992986R70627C- 
•»5215511P9?5057E“ 
-.70696553 56T0feqg- 
■.1175 6 X5839620*9- 
-.3302 901 7366229F- 
.1232'270783O531E^ 
- 36827653566036E- 


150751295058079. 


59525227153711E- 


. P860638291S156E-02 
.3575e06e795712E-01 
.51C63B55735C5eE-ul 
.31157577165291 E_^01 
- .32521537 596773E-01 
.25725026bb9195E-01 
Ifjb b7e76 3e2J7101E -01 
.75 9862130 18616E-02 
.177965271 'j9 U 6E-01 
.1 68B57617B2160E-01 
^ 1253070716ieB5E-(> l 
■. 21795 665862655E-02 
-. 36627653 56 ^33 E-01 
. 11388653681786 En,u 

-0 


-.13389798573057E-03 

.16271507517355E-03 

•.15e66302136599E-02 

.5052829o85299"5E-02 


. 23550319 
-.50775656 
-.95055056 
-.15009657 
-. 71 59652 5 
-.7025525 1 
-.1507512 9 
-.50535605 


1 5555562 

10171285 


359165 E- 
59E315E- 
527632E - 
51525 2E- 
6351606- 
9162 91E- 
595S'96t = 
■613li3T- 


^C^8 6£- 
3555351- 


COL 16 



-.103970628006926-01 




.602352077652<’3E-Ol 


3295O55365O552E-01 


175556625C 5032F-C1 


7C023036855657F-01 



.023051663*18626-02 


.210Q011*31B505E-01 


.25367506277302E-01 


.10256656150760E-C1 


.670520093P0236F-C2 


. 50525227155555E-OZ 


,k> I 


.10]71ZP530*677F-03 


.303C907F975055E-C1 
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TABLE 13. REAL TRANSFORMATION MATRIX (Y) 


sixsr 



6. 

07 

^ • 0 • 0 • 

d. 

0. 

0. ■ 

0. 

0, 




0. 

0. 

0. 0. 0. 

0. 

0, 

■d. 

6. 

0. 

0'.' 

J 


0. 

0. 

0 • 0 • 0 • 

0. 

0. 

0. 

0. 

0. 

0. 



0. 

0. 

^ • 0 • 0 • 

0. 

0. 

0. 

0. 

0. 

0. 



0. 

0. 

• 0 * Q *. 

0. 

0, 

0. 

0. 

0. 

0. 

' 


0. 

0 . 

0 * 0 • 0 » 

d". 

0. 

0. 

0. 

0. 

O'. 



0. 

0 . 

0. 0. 0. 

0. 

0. 

0. 

0. 

0. 

0". 



0. 

0. 

0 • Os 0 • 

0. 

D. 

0. 

0. 

0. 

0. 



0 . 

0. 

0. 0. 

0. 

0. 

'O'; 

0. 

0. 

O': 



0. 

0. 

0. 0. 0. 

d. 

D. 

0. 

0. 

0. 

0. 



^ • 

0 . 

0 • 0 • 0 • 

0. 

0. 

0. 

0. 

0. 

O'; 



0 . 

0. 

U • 0 • 0 • 

0. 

0. 

0. 

0. 

0. 

0. 



Q. 

0 . 

^ • 0 • 0 • 

0. 

0. 

"d. 

0. 

0. 

""0. 



0. 

0 . 

0 « 0 • 0 • 

0.' 

0. 

0". 

D. 

0. 

57 

' 


0 . 

0* 

0 • 0 t (J • 

"d. 

0. 

0. 

0. 

0. 

0. 


0. 

0 • 

0. 

Os Os Os 

0. 

O'. 

0. 

0. 

0. 

'0. 


=75n‘7T+'0'7 

.'lyFEy+or 

-.56676+06 

.r0'84E + 0'(r -.217FE+05 ."410‘22E+04" 

^75 4^6-6+ 0'3" 

.6"0 89TT+'05“ -,2T9'96 + 0'5 

■~7'7'9'50'E+04 

-.13216+04 


3'513'ET03'~ 

.24576+07 

-. 2304'6+D7 

. 151'9“E+07 
■■■■^r95'8't'+07" 

-.49676 + 06 .99726+'05 -.18446 + 05" 

72"4'9'6 =+d"4 

.T9T8E + 04" 

.3'8'2'9E+'0'4" 

■=".'13T2'6+0'5 

■,"6'06 3E + 04 

- . 

16105+04 

-.92 SITE + 06 

,16996+07 

.1155E+07 -.3749E+06 .69316 + 05' 

-79‘3S4'F + 04" 

-.76 80E + 0'4 

.13996+05 

.37376+04 

-.'70926+04 


6057E+04 . 


-.62046+06 .1290E+0r 

-.r'T24E+07 .■855'6E+06 -. 260'6'E+0'5' 

~T3'5'2'8 ?+ 0"5 '.T6'3'4 ET0'4‘" 

-T97586+'0'4 

.91356 + 04 =-. '25386+04 


'98 83 c + 04 

- . 4 5 2 2 r+ 0 5 

. 14196+06 

-.4771E+06 

.9749E+06 -.10616+07 . 577OE+Oi5“-.l‘2776+0'6 

"=737 39E + or 

'72 2 31E + 0'4" 

■=7 6'996E + 0'8 

.10596+05 

** s 

41796+04 

T^'BTE'E+O^' 

-.2463E+0y . 8277E+05 

-.27876+05 .54l‘6E + 06 -.5110“F+0'6" 

"~;r85T6+d'6" 

"■7'6'4"87E'+0’2" 

■=•738706 + 03" 

■~7r213'6 + 0'4‘ 

-.49556+04 


27076 + 04' 

-.ll‘*8E+04 

.3605F+04 

-.12126+05 

.4079E+'05' -.T295E+06 .20016 + o'5" 

-.10216+06 

-.9'49'5E + 01 

756656+02 

-,1776E'+03 

"'.72506+03 

- . 

1740F+O4 

.29096+07" 

.15576+06 

-.216'9E + 06 

.■4T33E+05"-.82 99E+04 .15346+04' 

"=';2077F+'0'3" 

■9t;-774'0'E+'06 

■■71'5 30E'+0"6" 

"";'3'l3 5 6 + 04 

-750406+03 


1340F+03 

-.173 2 r+D 7 

• ZO 6 tJ E + 06 

. 67'4'8'E +06 

-.42136+06 .84536+05 -.15636 + 05 7211'6'6+04' 

.'2'6"r2E + 0'6 

-.4'8T5 6''+0"6 

.2O71'E + 06 

. 5'170 6+'0"4" 

- * 

1355F"+04 ' 

. 9227E"+o6 

-.Q79 6'6 + 66~ 

.23636+06 

.51836+06 -.34846+06 .64366+05" 

-.87146+04 ,70336+04 

.27216+06 

-.45696+06 

.17466+05 


55336+08 

— i I i^09f:+0fc 

.53726 + 06" 

-.6"250E+06 

-.18006+06 .6593E+0b' -.3286E + 0'6 

.44496+05 

-,1413'6'+0'4 

,"8492 F + 04 

,21'82E+06 

-.3"777'E+05 


1471 = + '05 ■ 

.646'16+'05 

-,202 8E+"0"6 

.■68216 + 06 

-.99256 + 06' ,36986-^06 ,25536"+06" 

■^.152'46+0 6 

,534'2 6'+0'3 

-.31876+04 

■ .1002 6'+05 

,20916+05 


352‘26 + 06 

-.20336+65 

.63806+05 

-.21456+06 

.72225+06 -.18306+07 .38806+06 

.52326+06 

-.l"6816 + 03 

.10036+04 

-.31436+04 

.T2 8 5"E+05 


12206+05 

• 6^7 ^F + 03 

-. 27906 + 04 

.7361E+"04 

-.24786+05 .7868E+05 -.2508E+06" 

.19096+06 

.5769E+01 

-.344'1'E+62 

.10796+03 

-.44D5E+0‘f 


10506+08 

. 51TTET06 

-.16'0 56 + 07' 

.53966+07 

-.18176+08 .57706+06 .88376+07' 

-.52086+03 

,42296+04 

-.25236+05 

.79086+05 

-.32296+05 


77456+06 

0. 

0. 

0. 

Os Os Os 

"0"; 

'd'. ■ 

0. 

D. 

0. 

0. 


.4E)8 0E + O'5 

-.13156+05 

"". 45066 + 04" 

-.86186+03 .17306+03 -.31986+02' 

,'43296+01 

'-.225'4F'+d3 

,28056+02 

-.96666+02 

-.3986E+0'l 


81496+01 

-.1515E+05 

.r4'70E + 05' 

-.99706+04 

.33296+08 -.66836+03 , 12366+03" 

■^,16736+02 

■^72116'E + 0'3 

.62066+02 

-,3'533E + 02 

-.67396+02 


11006+01 

.^743Ft.04 

-78B5TSTGT 

. 10 436+ or 

-.63266+04 .21006+04 -.3e82E'+'03" 

.52566+02 

■^23256 + 02 

-.21276+03 

.32076+02 

-.52756+02 

-S 

55656+02 

-.009iTE + 0r 

.28556+04 -.60376+04 

.■68446+04' -.42266 + 04 ,13136 + 04 

-.17786+03 

-.5734F+02 

-.67606+01 

-.15226+03 

.19856+02 


31816+02 , 

#l787E*03 

-.56106+03 

.18866+04 

-.39176+04 .43756+04 -.2452F+04" 

.5577E + 0'3 

-.806"8'E+'02 

-.8878E+02 

-.2023'E+d2 

-. 14 5'9E+03 


33036+02 , 


I - . I e'TTE +‘CT2 .5‘8F6"e^"2 -. 1978E + C53 T66'6'0E+6’3 -.*r29'2EVo'4 .T3'07F+0V~-^5 30TE + 0 3 -.31 WeTo'2~-7278«E + 02''^.3'3^9E+02 -.30%9F+-3'2 -.75'3^E + 02 | 

^3 97 7E + DT_- .12V9F+0 2] . 4T9~e'E-t-0'2 T4T3E + 03_ .4 4 8 8 E>~ 0 3 - . 6 80 2 F ♦ 0 3 7 3'4 0 8 E ■*■ 0 3'~^ 1 9T3E » 0 2~^' 1 7 43 E ♦ 0 2'~~^ T 7 O' 1 FTo~2~ - 72 500E +0 2_- . 2 1 2 9 = ♦ CT2_| 

-.rr73'E+'04 -VQeE^E+'O'B .5478E+'0'3 -;i046£>03' '.2099E + 02 -73881E + 0r ;5255E + 00 . 1002 E+0 4’^Tl 854 E + 03“^790'4'E + 01 731258 + 00 -^7872E+0'0 ■ 

^v475yE+0‘3_ 7613 3E + 0 3 llT7E+'04_ .'541SE + 0 3 _2*108 7E+03 ' .■■2010E+02”^”72 722 £>0 1“ ^72 1 90E y03“~V4 3 1 3E + 03 ~7l 735'F+03_f76 5 2 9E + 0l'_ . n7<i'F+01 _ 
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